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Abstract 

In this paper, we study feedback optimization problems that maximize the users' signal to interference plus 
noise ratio (SINR) in a two-cell MIMO broadcast channel. Assuming the users learn their direct and interfering 
channels perfectly, they can feed back this information to the base stations (BSs) over the uplink channels. The 
BSs then use the channel information to design their transmission scheme. Two types of feedback are considered: 
analog and digital. In the analog feedback case, the users send their unquantized and uncoded CSI over the uplink 
channels. In this context, given a user's fixed transmit power, we investigate how he/she should optimally allocate 
it to feed back the direct and interfering (or cross) CSI for two types of base station cooperation schemes, namely, 
Multi-Cell Processing (MCP) and Coordinated Beamforming (CBf). In the digital feedback case, the direct and 
cross link channel vectors of each user are quantized separately, each using RVQ, with different size codebooks. 
The users then send the index of the quantization vector in the corresponding codebook to the BSs. Similar to the 
feedback optimization problem in the analog feedback, we investigate the optimal bit partitioning for the direct and 
interfering link for both types of cooperation. 

We focus on regularized channel inversion precoding structures and perform our analysis in the large system 
limit in which the number of users per cell (K) and the number of antennas per BS (N) tend to infinity with 
their ratio 13 = ^ held fixed. We show that for both types of cooperation, for some values of interfering channel 
gain, usually at low values, no cooperation between the base stations is preferred: This is because, for these values 
of cross channel gain, the channel estimates for the cross link are not accurate enough for their knowledge to 
contribute to improving the SINR and there is no benefit in doing base station cooperation under that condition. We 
also show that for the MCP scheme, unlike in the perfect CSI case, the SINR improves only when the interfering 
channel gain is above a certain threshold. 
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I. Introduction 

A. Background 

Recently, many applications that require high data rates such as high quality video streaming and huge volume 
data transfers through wireless communication systems have emerged. MIMO communication systems have arisen 
as a promising candidate to support this requirement and have been adopted for existing and future wireless 
communication standards such as in IEEE 802.1 In and 4G networks. Current MIMO technological advancements 
can be considered as the results of research works started about fifteen years ago. So far, there has been a 
considerable amout of work focusing on single user and single-cell multiuser MIMO systems. Only recently, 
researchers have started to put more attention to investigate how to maximize data rates in multi-cell MIMO 
networks, particularly in the downlink [1, and references therein]. 

The main challenge that limits the spectral efficiency in the downlink of multi-cell networks, besides intra-cell 
interference, is the inter-cell interference (ICI). The conventional approach to mitigate this interference is to use 
spatial reuse of resources such as frequency and time [1]. The move towards aggressive frequency or time reuse 
will cause the networks to be interference limited especially for the users at the cell edge. The current view is to 
mitigate ICI through base station (BS) cooperations. Within this scheme, the BSs share the control signal, channel 
state information (CSI) and data symbols for all users via a central processing unit or wired backhaul links [2]. 

It has been established in [3]-[8], to name a few, that MIMO cooperation schemes provide a significant increase 
in spectral efficiency compared to conventional cellular networks. BS cooperation can be implemented at different 
levels [1]. In the Multi-Cell Processing setup, also known as Network MIMO or Coordinated Multi-Point (CoMP) 
transmission, the BSs fully cooperate and share both the channel state information (CSI) and transmission data. 
This full cooperation requires high capacity backhaul links which are sometimes not viable in practical settings. 
To alleviate this requirement, only CSI (including direct and interfering channels) is shared amongst base stations 
in the interference coordination scheme [1]. Several works have addressed coordinated beamforming and power 
control schemes to improve the spectral efficiency in interference-limited downlink multi-cell networks. Detailed 
discussions regarding these topics can be found in [1] and references therein. 

In both base station cooperation schemes, the CSI at the base stations plays an important role in maximizing the 
system performance. The base stations use this information to adapt their transmission strategies to the channel 
conditions. The benefit of having CSI at the transmitter (CSIT) with respect to the capacity in single and multi- 
cell multi-antenna systems is nicely summarized in [9], [10]. However, these advantages are also accompanied by 
the overhead cost for the CSI acquisition via channel training and feedback in frequency division duplex (FDD) 
systems. It needs to scale proportionally to the number of transmit and receive antennas and the number of users 
in the system in order to maintain a constant gap of the sum-rate with respect to the full CSI case [11]. Moreover, 
in practical systems, the backhaul-link capacity for CSI and user data exchanges and feedback-link bandwidth are 
limited [2]. Considering the CSI signaling overhead from channel training and CSI feedback, references [12], [13] 
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(see also [14]) suggested that the conventional single-cell processing (SCP) without coordination may outperform 
the cooperative systems, even the MCP scheme. In this paper, to reduce the complexity in the analysis, we ignore 
the (important) constraints of limited backhaul-link and CSI training overhead. We assume a perfect CSI training 
so that all users know their CSI perfectly. We focus on studying how to allocate feedback resources, that depend 
on the feedback schemes, to send the CSI for the direct channel and interfering (cross) channel to BSs so that 
the users' SINR are maximized. Two feedback schemes are considered in our study: the analog feedback scheme, 
introduced in [15] and the limited (quantized) feedback via random vector quantization (RVQ), introduced in [16]. 
In the analog feedback scheme, each user sends its unquantized and uncoded channel state information through 
the uplink channel. Hence, we ask the question, for a given uplink power constraint, what fraction of this uplink 
power is allocated optimally to transmit the direct and interfering channel information? For the digital feedback 
scheme, the number of feedback bits determines the quality of the CSI. Hence, we can ask, how many bits are 
optimally needed to feedback the direct and cross CSI? 

B. Contributions 

The main goal of this paper is to optimize and investigate the effect of feedback for MCP and CBf cooperation 
schemes under analog and quantized feedback (via RVQ). We consider a symmetric two-cell Multi-Input Single- 
Output (MISO) network where the base stations have multiple antennas and each user has a single antenna. We 
assume that the users in each cell know their own channel perfectly: they feed back this information through the 
uplink channel and the base stations form the users' channel estimates. The BSs use these estimates to construct a 
regularized channel inversion (RCI) type beamformer, also called regularized zero-forcing (RZF), to precode the data 
symbols of the users. The precoders follow the structures proposed in [17]. Unlike [2], [14], we assume several users 
are simultaneously active in each cell so that the users experience both intra- and inter-cell interference. To mitigate 
ICI through base station cooperation, we consider both fuU cooperation (MCP) and interference coordination via 
CBf. 

Our contributions can be summarized as follows. First, under both feedback models and both cooperation 
schemes, we derive the SINR expression in the large system limit, also called the limiting SINR, where the number 
of antennas at base stations and the number of users in each cell go to infinity with their ratio kept fixed: As our 
numerical results will show, this is indicative of the average performance for even finite numbers of antennas. 
Then, we formulate a joint optimization problem that performs the feedback optimization for both feedback models 
and both cooperation schemes and finds the optimal regularization parameter of the corresponding RCI-structured 
precoder. The regularization parameter is an important design parameter for the precoder because it controls the 
amount of interference introduced to the users. Optimizing this parameter, as discussed later, will allow the precoder 
to adapt to the changes of the CSIT quality and consequently produces a 'robust beamformer' . 

We analyze the behavior of the maximum limiting SINR as a function of the cross channel gains and the 
available feedback resources, and identify, for both the analog and quantized feedback models, regions where SCP 
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processing is optimal. We also show that whereas in the perfect CSI case, MCP performance always improves 
with epsilon, this only occurs after a certain threshold is crossed in both analog and limited feedback cases. 
Parts of this work appeared in [18], [19], but without the proofs. 

C. Related Works 

In the last decade, there has been a large volume of research discussing feedback schemes in multi-antenna 
systems. A summary of digital feedback (also known as limited or finite-rate feedback) schemes in multi-antenna 
(also single-antenna) and multi-user systems in the single-cell setup can be found in [11]. Since the optimal 
codebook for the limited feedback is not known yet [16], [20], [21], the use of RVQ, which is based on a random 
codebook, as the feedback scheme becomes popular Furthermore, the RVQ-based system performance analysis is 
also more tractable. In multiple-antenna and multiuser systems, works on the analog feedback commonly refer to 
[15] (sometimes [22]). 

The paper by Jindal [20] sparked the use of RVQ in analyzing broadcast channels. Considering a MISO broadcast 
channel with a zero-forcing (ZF) precoder and assuming that each user knows its own channel, the main result in 
the paper is that the feedback rate should be increased linearly with the signal-to-noise ratio (SNR) to maintain the 
full multiplexing gain. Caire et al. in [21] investigate achievable ergodic sum rates of BC with ZF precoder under 
several practical scenarios. The CSI acquisition involves four steps; downlink training, CSI feedback, beamformer 
selection and dedicated training where each user will try to estimate the coupling between its channel and the 
beamforming vectors. They derive and compare the lower bound and upper bound of the achievable ergodic sum- 
rate of the analog feedback as in [15] and RVQ-based digital feedback under different considerations, e.g., feedback 
transmission over AWGN and MAC channel, feedback delay and feedback errors for the digital feedback scheme. 
A subsequent work by Kobayashi et al. in [23] studies training and feedback optimizations for the same system 
setup as in [21] except without dedicated training. The optimal period for the training and feedback that minimized 
achievable rate gap (with and without perfect CSI) are derived under different scenarios as in [21]. The authors 
also show that the digital feedback can give a significant advantage over the analog feedback. In the same spirit as 
[20] , reference [24] discusses the feedback scaling (as SNR increases) in order to maintain a constant rate gap for 
a broadcast channel with RCI precoder. The analysis has been done in the large system limit since the analysis the 
finite-size turns out to be difficult [20]. Moreover, besides analyzing for the case K = N, as in [20], the authors 
also investigate the case K < N. 

While channel state feedback in the single-cell system has received a considerable amount of attention so far, 
fewer works have addressed this problem in multi-cell settings. The effect of channel uncertainty, specifically 
the channel estimation error, in the multi-cell setup is studied in [25], [26]. In [26], the authors conclude that 
when channel estimates at one base station contain interferences from the users in other cells, also called as pilot 
contamination phenomenon, the inter-cell interference increases. Thus, this phenomenon could severely impact 
the performance of the systems. Huh et al. in [25] investigate optimal user scheduUng strategies to reduce the 
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feedback and also the effects of channel estimation error on the ergodic sum-rate of the clustered Network MIMO 
systems. They consider the ZF precoder at the base stations and derive the optimal power allocation that maximizes 
the weighted sum-rate. In deriving the results, it is assumed that the BSs received perfectly (error-free) the CSI 
fed back by the users. The overhead caused by the channel training is also investigated and they observe that 
there is a trade-off between the number of cooperating antennas and the cost of estimating the channel. Based on 
the trade-off, the optimal cooperation cluster size can be determined. By incorporating the channel training cost, 
no-coordination amongst the base stations could be preferable. The same conclusion is also obtained in [12], [13]. 

For the interference coordination scheme, [2] presents the RVQ-based limited feedback in an infinite Wyner 
cellular model using generalized eigenvector beamforming at the base stations. The work adopts the intra-cell 
TDMA mechanism where a single user is active in each cell per time slot. Each user in each cell is also assumed 
to know its downlink channel perfectly. Based on that system model, an optimal bit partitioning strategy for direct 
and interfering channels that minimizes the sum-rate gap is proposed. Explicitly, it is a function of the received 
SNR from the direct and cross links. It is observed that as the received SNR from the cross link increases, more bits 
are allocated to quantize the cross channel. A better quality of the cross channel estimate will help to reduce the 
inter-cell interference. The authors also show that the proposed bit partitioning scheme reduces the average sum- 
rate loss. Also in the interference coordination setting, [14] takes into account both CSI training and feedback in 
analyzing the system what they called the inter-cell interference cancellation (ICIC) scheme. In ICIC, the precoding 
vector of a user is the projection of its channel in the null-space of the others users' channels in other cells so 
that the transmission from this user will not cause interference to the users in other cells. The work also assumes 
the intra-cell TDMA and presents the training optimization and feedback optimization for both analog and digital 
feedback (RVQ). Based on that system setup, the most interesting result is that the training optimization is more 
important than the feedback optimization for the analog feedback while the opposite holds for the digital feedback. 

For different levels of cooperation, i.e., MCP, CBf and SCP, [17] investigates an optimization problem to minimize 
the total downlink transmit power while satisfying a specified SINR target. The authors derived the optimal transmit 
power, beamforming vectors, cell loading and achieved SINR for those different cooperation schemes in a symmetric 
two-cell network. The resulted optimal beamforming vectors have a structure related to RCI. 

The current work is closely related to [17] in the sense we use the same cooperative schemes and precoder 
structure. We extend the work by analyzing the optimal feedback strategies for analog and digital feedback under 
MCP and CBf schemes. The results in this work are obtained by performing the analysis in the large system limit 
where the dimensions of the system i.e., the numbers of users and transmit antennas tend to infinity with their 
ratio being fixed. The large system analysis mainly exploits the eigenvalue distribution of large random matrices. 
For examples, it has been used to derive the asymptotic performance of linear multiuser receivers in CDMA 
communications in early 2000 (see [27]), single-cell broadcast channels with RCI for various channel conditions 
[24], [28]-[30], base station cooperations in downlink multi-cell networks (see e.g., [17], [25]). The asymptotic 
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performance measure becomes a deterministic quantity and can have close-form/compact expressions. Hence, it 
can be used to derive the optimal parameters for the system design. Moreover, it can provide a good approximation 
of, hence insights on, the performance of the finite-size (or even small-size) systems. 

Similar to [2] and [14], we perform the feedback optimization in interference-coordination scheme (CBf). As 
in [14], we also investigate the feedback optimization for the analog and digital feedback schemes. However, 
different from those works, we do not assume the intra-cell TDMA in each cell, and hence each user experiences 
both intra-cell and inter-cell interference. We also consider a different type of precoder i.e., the RCI. Moreover, we 
also analyze the feedback optimization for different level of cooperations between the base stations, including the 
MCP setup, and try to capture how we allocate resources available at the user side as the the interfering channel 
gain varies. 

D. Paper Organization and Notation 

The rest of the paper is structured as follows. The system model is described in Section II. It starts with the 
channel model, and the expressions of the transmit signal, precoder and the corresponding SINR for each MCP and 
CBf. In the end of the section, the feedback schemes and true channel model in term of the channel estimate at the 
BSs and the channel uncertainty for the analog and digital feedback are presented. The main results for the noisy 
analog feedback and digital feedback and for different types of coordination are discussed in Section III and IV, 
respectively. In each section, we begin by discussing the large system result of the SINR for the MCP and CBf and 
then followed by deriving the corresponding optimal feedback allocation; optimal (uplink) power for the analog 
feedback and optimal bit partitioning for the digital feedback. The optimal regularization parameter for the RCI 
precoder is also derived for both types of feedback and cooperation. The end of each section provides numerical 
results that depict how the optimal feedback allocation and the SINR of each user behave as the interfering channel 
gain varies. In Section V, we provide some numerical simulations that compare the performance of the system 
under the analog feedback and digital feedback. The conclusion are drawn in the Section VI and some of the 
proofs go to the appendices. 

Throughout the paper, the following notations are used. E[-] denotes the statistical expectation. The almost sure 
convergence, convergence in probability, and mean-square convergence are denoted by respectively. 
The partial derivative of / with respect to (w.rt.) x is denoted by The circularly symmetric complex Gaussian 
(CSCG) vector with mean // and covariance matrix S is denoted by CM {/J,, XI). \a\ and 3? [a] denote the magnitude 
and the real part of the complex variable a, respectively. || • || represents the Euclidean norm and Tr (•) denotes 
the trace of a matrix. I^v and On denote an x identity matrix and a 1 x zero entries vector, respectively. 
(•)-^ and (•)^ refer to the transpose and Hermitian transpose, respectively. The angle between vector x and y is 
denoted by Z(x, y). LHS and RHS refer to the left-hand and right-hand side of an equation, respectively. 
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II. System Model 

We consider a symmetric two-cell broadcast channel, as shown in Figure 1, where each cell has K single 
antenna users and a base station equipped with antennas. The channel between user k in cell j and the BS in 
cell i is denoted by row vector h^j j where hj^j j ~ CM{Q,1n) and h^jj ~ CAA(0,elAr), for j = 1,2 and j = 
mod {j, 2) + 1. We refer to the .,■ j as direct channels and h;. as cross or "interfering" channels. We find it 
useful to group these into a single channel vector = [h^ j i j 2]- 

We consider an FDD system and assume that the users have perfect knowledge of their downlink channels, 
hfc j and ^j. Each user feeds back the channel information to the direct BS and neighboring BS through the 
corresponding uplink channels. The BSs estimate or recover these channel states and use them to construct the 
precoder. 

The received signal of user k in cell j can be written as 

yk,j = hfc,j,ixi + hfcj^2X2 + nfcj, 

where Xj G C^^^, z = 1,2 is the transmitted data from BS i, and ~ CA/'(0,cj^) is the noise at the user's 
receiver. The transmitted data Xj depends on the level of cooperation assumed, and will be described in more details 
in Sections II-A and II-B: we restrict ourselves to linear precoding schemes, more specifically RCI precoder. We 
assume each BS's transmission is subject to a power constraint E [||xj|p] = Pj. In the MCP case, we relax this 
constraint to a sum power constraint so that E [||x|p] = Y^^=i — ^t- In the analysis, we assume Pi = P2 = P 
and denote 7^ = P/o-"^- 

As already mentioned, in practical scenarios, perfect CSI is difficult to obtain and the CSI at the BSs is obtained 
through feedback from the users. We are particularly interested in the channel model where we can express the 
downlink channel between the user k in cell j and BS i as 

hfcj,i = -x/ 4>k,j,i^k,j,i + hfcj,j, (1) 

where j j represents the channel estimate, and j j the channel uncertainty or estimation error. The channel 
estimates are used by the BSs to construct the precoder. 

The transmitted signal, precoder and SINR for each user for each cooperation scheme will be presented in the 
following subsections. 

A. MCP 

As previously mentioned, in the MCP, both BSs share the channel information and data symbols for all users in 
the network. Therefore, we may consider the network as a broadcast channel with 2A^ transmit antennas and 2K 
single antenna users. The BSs construct the precoding matrix by using the users' channel estimates. In this work, 
we consider the RCI precoding, for which the precoding or beamforming vector for user k in cell j, w^j, can be 
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Fig. 1. System model 



written as [31] 



where h^j = [h^ j i hfc _,- 2] and H = [h^^ h^^^^ • • • ^ hf^2 ^22 ' ' ' ^K2\^- "^^^ transmitted data vector can 
be expressed as 

2 K 

j=l k=l 

where s^j ~ CJ\f{0, 1) denotes the symbol to be transmitted to user k in cell j. It is also assumed that the data 
symbols across the users are independent, i.e, E[ss^] = I2N, with s = [si 82]^ and Sj = [sij S2j ■ ■ ■ sxj]'^- c is 
a scaling factor ensuring the total power constraint is met with equality: 

Pt 



The received signal at user k in cell j can be written as 



Vkj = hfc ,jX + rifc J = chfc J ( H-^H + al27v ) H^" s + nj. 



where follows the channel model (1) with j = [h^ j 1 hfcj 2]- The term j and j are obtained from 
H and s by removing the row corresponding to user k in cell j respectively. Hence, the SINR for user k in cell 
j can be expressed as 

2 



-1 



SINR 



hfc,,- (H^H + al2jv) h^l, 



c^hfc,,- (H^H + al2N) H|^,Hfc,i ( H^H + alsiv j h^^, + <yl 



(2) 
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B. Coordinated Beamforming 

In this scheme, the base stations only share the channel information, so that, for cell j, Xj can be expressed as 

K 
k=l 

where as in the MCP case s^j ~ CJ\f{0, 1) denotes the symbol to be transmitted to user k in cell j. The constant 
Cj is chosen to satisfy the per-BS power constraint, that is, E [||xj|p] = Pj. Hence, c| = ^^'^^^^^2 - We let 

( - V- 

^kj = 1 OH-N + ^ hf^^ jhi^rnj j ^k,j,j, 

\ {hm)^(k,3) j 

which is an extension of regularized zero-forcing to the coordinated beamforming setup [17]. Note that designing 
the precoding matrix at BS j requires local CSI only (the h,fc j j from BS j to all users, but not the channels from 
the other BS to the users). The SINR of user k in cell j can be expressed as 

2^ (^r\^k,j,r^k'r\ +(yd 

where, once again, j and h^ j j/ follow (1). 

C. Analog Feedback through AWGN Channel 

In the analog feedback scheme, proposed in [15], each user feeds back the CSI to the base stations using the linear 
analog modulation. Since we skip quantizing and coding the channel information, we can convey this information 
very rapidly [15]. We also consider a simple uplink channel model, an AWGN channel. A more realistic multiple 
access (MAC) uplink channel model could be a subject for future investigation. Each user in cell j feeds back its 
CSI hfc orthogonally (in time). Since each user has to transmit 2N symbols (its channel coefficients), it needs 
2kN channel uses to feed back the CSI, where k > 1. User k in cell j sends 

hfcjA|, (4) 

where Aj is a diagonal matrix such that the first diagonal entries are equal to Aji and the remaining diagonal 
entries are equal to Aj2, with \jj = IukPu, Xjj = 2e~^(l — v)kPu and Pu is the user's average transmit power per 
channel use. Equation (4) satisfies the uplink power constraint E[||h,fc ^ A| |p] = 2kNPu. Thus, the power allocated 
to feedback the direct and interfering channel is controlled by G [0, 1]. We should note that in (4), it is assumed 
that K is an integer. If kN is an integer, we can modulate the signal (4) with 2N x 2kN spreading matrix [15], 
[21] and the analysis presented below still holds. 

Now, let be, i = 1,2, ■ ■ ■ , 2N, be the Ah element of hkj, Xe be the corresponding element on the diagonal of 
A, and ee = E[6^6|]. When this channel coefficient is transmitted, the signal received by the coordinating BSs is 
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In 



where G C^^^^ ~ CN{Q,all) is the noise vector at the coordinating BSs and l^r is a column vector of length 
N with all 1 entries. Using the fact that the path-gain from the users in cell j to BS j is e, the MMSE estimate 
of becomes 

hi = \/%eip^ [A^e^pp^ + alhN] ^ Ye, 

and its MMSE is cr^^ = — A^e^p-^ [A^e^pp-^ + (T^l2Ar] ^ P- We should note that {b^} are mutually independent. 
By using the property of MMSE estimation, we can express j j as 



(5) 



where j j represents the channel estimate, and j the channel uncertainty or estimation error. Note that the 
entries of each vector j ^ and j j are independent and identically distributed (i.i.d) and distributed according 

to CM{0,ujji) and CM{0,5ji), respectively, where 



l+(l-i/)7„ ' 



J = i 



e(l-i/)7„ 
l+(l-i/)7„ 



J = I 



(6) 



and 7„ = 27uk(1 + e) with 7^ = NPu/a"^. The channel estimates are used by the BSs to construct the precoder. 
Since each 6ij and ujij are identical for all users then we denote 6d = 6c = ^fj^^d = ^jj and ujc = ujfj. From 



(6), it follows that LOd = I — dd and ujc 



6r. 



D. Quantized Feedback via RVQ 

In the digital feedback case, user k in cell j uses B^jj and B/^ jj bits to quantize or feedback the direct and 
interfering channels, respectively. The total number of feedback bits is assumed to be fixed. It is also assumed 
that each user has different codebooks: Ukjj with size 2^'= and U/^jj with size 2^'' J 3, to quantize the direct 
and interfering channel, respectively. Moreover, these codebooks are different for each user. In this work, Bf^ j j 
is the same for all users and B^jj = Bd,yk,j = 1,2. Similarly, Bf. jj = Bc,'ik,j = 1,2. The total number of 
feedback bits is denoted by Bt, where Bt = Bd + Be- 

Since the optimal codebook design for the quantized feedback is not known yet, therefore in this paper, for 
analytical tractability, we consider the well known RVQ scheme. As suggested by its name, RVQ uses a random 
vector quantization codebook where the quantization vectors in the codebook are independently chosen from the 
isotropic distribution on the A^-dimensional unit sphere [16], [20]. The codebook is known by the base station 
and the user. The user quantizes its channel by finding the quantization vector in the codebook which is closest 
to its channel vector and feedbacks the index of the quantization vector to the BSs. We should note that only 
the channel direction is quantized. Most of the works that employ RVQ for the feedback model assume that only 
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channel direction information is sent to the BSs. As mentioned in [20], the channel norm information can also be 
used for some problems that need channel quality information (CQI) such as power allocation across the channel 
and users scheduling [32]. 

The user k in cell j finds its quantization vector for the channel j j according to 



= arg max 



\^k,j,i'^k,j,i\ 



The quantization error or distortion r| j ■ is defined as 

||hfc -Ufe -ll^ 

It is a random variable whose distribution is equivalent to the minimum of 2^*=^ ' beta random variables with 
parameters — 1 and 1 (see [20], [33]). Each realization of Tkj^i is different for each user even though the users 
have the same amount of feedback bits. 

Having obtained Uk,j,i, each user then sends its index in the codebook and also the channel magnitude ||hfe j|| 
(see also [32]). By assuming that the BSs can receive the information perfectly, the channel estimate at the BS 
can be written as 

— llufc j'^j. (7) 

Note that h^j^i has the same statistical distribution as j j i.e., h^j ^ ~ CAA(0, ejilN), where eji = 1 when i = j 
and otherwise, eji = e. 

From [20], [34], we can model j j as follows 



^k,j,i — yl ''"fcj-^jhfcj-^j + Tfc j^j||h/jj^j||z/jj^j, (8) 

where z^j^i is isotropically distributed in the null-space of Uk,j,i and is independent of j j. Moreover, z^j j can 
be rewritten as 

vi. ,-11- 

Zk,j,i 



i|vfcj,ini II 

where Ilr is the projection matrix in the column space of j, 11- = In—-w^'^—^^ and v^j j ~ CM{0,1n) 
is independent of hkj.i- It is clear that the channel model (8) has the same structure as (1) with (l)k,j,i = 1 — 7"| ^ ^ 
and h = Tkj,i\\hkj^i\\zkj^i. 

E. Achievable and limiting sum-rate 

Besides SINR^ j, another relevant performance measure is the achievable rate. For the user k at cell j, it is 
defined as 

=log2(l + SINRfcj), (9) 
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It is obtained by treating the interferences as noise or equivalently performing single-user decoding at the receiver. 
Observing (9), it is obvious that there is a one-to-one continuous mapping between the SINR and the achievable 
rate (see also [35]). The total sum-rate, or just the sum-rate, can then be defined as follows 

2 K 

Rsum = ^^Rkj. (10) 

j=l k=l 

As shown later in Section III and IV, as K,N ^ oo, we have 

SINRfcj - SINR°° ^ 0, (11) 

where SINR°° is a deterministic quantity and also called the limiting SINR. It is also shown that the limiting 
SINR is the same for all users. By using the result (11) and based on the continuous mapping theorem [36], the 
following 

holds (see also [24]) where the limiting achievable sum-rate can be expressed as R^^ = /31og2(l + SINR°°). For 
the numerical simulations, we also introduce the normalized sum-rate difference, defined as 

Aiisum - — x¥Ti^n — ' ^^^^ 

2N L^sumJ 

that quantifies the sum-rate difference, [Rsum] — R'^m' compared to the (actual) finite-size system average 
sum-rate. 

III. MCP AND CBf with Noisy Analog Feedback 

In this section, we will discuss the large system results and feedback optimization for the MCP and CBf by using 
the analog feedback model discussed in Section II-C. First, the large system limit expression for the SINR is derived. 
Then, the corresponding optimal regularization parameter that maximizes the limiting SINR is investigated. Finally, 
the optimal u that maximizes the limiting SINR that already incorporates the optimal regularization parameter will 
be discussed. 

A. MCP 

We start with the theorem that states the large system limit of the SINR (2). 



Theorem 1. Let p^^p = (Ufi + ujc) ^a/N and g{(3,p) be the solution of g{j3,p) = yp + j+^^^j ■ large 
system limit, the SINR of MCP given in (2) converges in probability to a deterministic quantity given by 

l + ^(l + ff(/3,p^,,, ))2 

7e + (1 + 9{/3, Pm.af)) 



SINR^c,AF = 7e5(/3, Pmaf) ^- (13) 
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where the effective SNR 7e is expressed as 

7e = 7 ; r = ^ ; (^4) 

^d + ^c + Y^ Sd + ^c + Y, 

Proof: See Appendix II- A ■ 
It is obvious from above that the limiting SINR is the same for all users in both cells. This is due to the channel 
statistics of all users in both cells are the same. The channel uncertainty, captured by uj, and 5,, affects the system 
performance (limiting SINR) via the effective SNR and regularization parameter /9m,af- 

As discussed previously, the (effective) regularization parameter Pm,af controls the amount of interference intro- 
duced to the users and provides the trade-off between suppressing the inter-user interference and increasing desired 
signal energy. The optimal choice of /3m,af that maximizes (13) is given in the following. 

Corollary 1. The optimal Pm.af that maximizes SINR'^^^.p^^p is 

P*M.AF = — ' (15) 

7e 



and the corresponding limiting SINR is 



SINRl'cZ, = gi/3,pl,,). (16) 



Proof: The proof follows easily from [37]. ■ 
It is interesting to see that the limiting SINR expression with p* becomes simpler and it depends only the 
cell-loading (/3) and the effective SNR. Clearly from (14), 7e is a function of the total MSE, 6t = 6d + 5c, that can 
be thought as a reasonable measure of the CSIT quality. Thus, adjusts its value as 5t changes. Now, from 
(14), it is obvious that 7e is a decreasing function of 5t. As a result, pl^,^^ is increasing with 5t. In other words, if 
the total quality of CSIT improves then the regularization parameter becomes smaller. In the perfect CSIT case, 
i.e., 5t = 0, and in the high SNR regime, p* ^p goes to zero and we have the ZF precoder. 

Now, we will investigate how to allocate v to maximize the limiting SINR (16), or equivalently g{l3, pl^/^^). z/ is 
captured by 7e (or p^af) ^d- It can be shown that g is decreasing (increasing) in Pmaf (7e)- Then, for a fixed 
/3 the limiting SINR is maximized by solving the following optimization problem 

e- 5c + l- 5d 
max 7e = 

As mentioned earlier, 7e is a decreasing function of 6t. Thus, the optimization problem above can be rewritten as 

mill 6t = 6d + dc = + — — • (17) 

From the above, it is very interesting to note that the optimal u that maximizes SINR^^ is the same as the one 
that minimizes the total MSE, 5t. 

It is easy to check that the optimization problem above is a convex program and the optimal u, denoted by v*, 
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can be expressed as follows 

0, V~e>% + 1 



1, ^/^<^ (IS) 

l+^(l-v/g) 
1+v^ ' 



7„ + l 

Otherwise. 



As a result, for -y/e < - , the BSs should not waste resources trying to learn about the "interfering" channel 
states. In this situation, the coordination breaks down and the base stations perform SCR The completely opposite 
scenario, in which the BSs should not learn the "direct" channels, occurs when y/e > 7^ + 1. Clearly, this can 
only happen if e > 1. When ^/e > 7^ + 1, the BSs also perform SCP but each BS transmits to the users in the 
neighboring cell. 

We end this subsection by characterizing the behavior of 7e (equivalently SINR^ct ), after optimal feedback 
power allocation, as the cross channel gain e varies. This also implicitly shows how the total MSE, 5t, affects the 
limiting SINR. Let 7^ = jj^- We analyze the different cases in (18) separately. 

1) -y/e < ^-qiY- "^^^^ is ^^^^ when the BSs perform SCP for the users in their own cell. For fixed ^u, this 
inequality is equivalent to e < e^^x' where e^^x ^ satisfies y^e|^ = ^ (^_)_^scp ■ Now, by taking the first 
derivative ^ and setting it to zero, the (unique) stationary point is given by 

SCP ^ ^ _ I 

If £ [Oi \/ ^max] ' it is easy to check that the limiting SINR is increasing until e = e^" and then decreasing. 

If yjldlu > 1 then e^p^ < 0, or equivalently, ^ < 0. Consequently, for this case, the limiting SINR is decreasing 



in e. Moreover, > yefnax if following condition holds 

V Idlui'i - 27d - lu) > (27d7« - 7(i - 7«)> (19) 

in which case ^ > 0, which implies that the limiting SINR always increases over e for this case. 

This behavior of 7e as a function of e can be intuitively explained as follows. When u = 1, the total MSE 
is = [i+e)^ +1 ~^ ^' where the first and second terms are 6^ and Sc, respectively. As e increases, 64 decreases 
whereas 6c increases. This shows that there is a trade-off between the quality of the direct channel and the strength 
of the interference. The trade-off is also influenced by parameters 7^^ and As shown in the analysis, when 
y/'jdju > 1> the effect of cross channel to the limiting SINR dominates. In contrast, if the condition in (19) is 
satisfied, the effect of the quality of the direct channel {6t) becomes dominant. If the aforementioned conditions 
do not hold, 6t causes the SINR to increase until e^p"" and after that the interference from the cross channel takes 
over as the dominant factor, thereby reducing the limiting SINR. 

2) 7« + l > \/e > Y^-' Here, the BSs perform MCR By taking ^ in that interval of e, it can be shown that we 
have a unique stationary which we denote as y^e^. We can also show that 7e is a convex function for e G [0,1] and 
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is increasing for e > 1. Thus, if < ^/€^ < 7m + 1, the limiting SINR will decrease for ^/e G a/^S 
and increase after that; Otherwise, the limiting SINR increases in the region. Here, for ^/e £ [^-^,1], we still 
can see the effect of the trade-off within 6t to the limiting SINR as e changes. In that interval, the quality of 
the direct channel becomes better as e increases; However, that of the cross channel decreases and this affects 
the SINR badly until e"p. After this point, the improvement in the quality of the direct channel wiU outweigh the 
deterioration of that of the cross channel, causing the SINR to increase. 

3) ^fi > 7m + 1." In this case, each BS performs SCP, but serves the other cell's users. We can establish that 
^ > 0. Hence, for this case, the limiting SINR is increasing in e. 

B. Coordinated Beamforming 

Theorem 2. Let p^af = jj, cind let Ta be the solution of the following cubic equation 

Ta = w-^ • (20) 

In the large system limit, the SINR of the coordinated beamforming given in (3) converges almost surely to a 
deterministic quantity given by 



fTA 



PcAF + (i+a;,r^)2 ^ (i+uar^y^ 



(21) 



+^d + Sc+ (i+^,r,4)2 + (i+tXr^)2 
Proof: See Appendix III- A ■ 
Similar to the MCP case, the limiting SINR expression (21) is the same for all users. By comparing (15) and 
(22), it is also interesting to see that Pcaf — Pm,af a given a. The optimal /Jc,af that maximizes the limiting 
SINR (21) is given in the following. 

Corollary 2. The limiting SINR (21) is maximized by choosing the regularization parameter according to 

PcAF = I3(-+Sd + S,] . (22) 



Id 



and the corresponding limiting SINR is 



sinr^Taf = ^dr*A, (23) 



where T\ is Ta with Pqaf = pIaf- 

Proof: Let 7e// = /3 (7^ ^ + Sd + 4) and ^ = ^^/j^^^y. + (i+^r^- It is easy to show that 

aSINRgg,,p 7e//-/>c.AF 

OPcAF [jeff + OPcAF 

where = -2/3 ( „ + > with < is given by (61). Thus, it follows that 

OpcAF '^OPCAF \(l+l^drA)^ J Opc.AF ° J \ ^ ' 

Pc.AF = 7e// is the unique stationary point and the global optimizer. Plugging back Pc af into (21) yields (23). ■ 
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Similar to the MCP case, the corollary above shows that the optimal regularization parameter adapts to the 
changes of CSIT quality and it is a decreasing function of 6t. 

Finding u that maximizes the limiting SINR of the CBf is more complicated than in the MCP case. It is 
equivalent to maximizing w^F, such that u G [0, 1]: this is a non-convex program. However, the maximizer u* is 
one of foUowings: the boundaries of the feasible set (z/ = {0, 1}) or the stationary point, denoted by which is 
the solution of 

"° - ^ (25) 



The point z/ = can be eliminated from the feasible set since the derivative of the limiting SINR with respect to 
u at this point is always positive. 

C. Numerical Results 

Since propagation channels fluctuate, the SINR expressions in (2) and (3) are random quantities. Consequently, 
the average sum-rates are also random. Figure 2 illustrates how the random average sum-rates approach the limiting 
sum-rates as the dimensions of the system increase. This is quantified by the normalized sum-rate difference which 
is defined in (12). The average sum-rate is obtained by averaging the sum-rates over 1000 channel realizations. The 
optimal regularization parameter and power splitting obtained in the large system analysis are used in computing 
the limiting and average sum-rates. We can see that as the system size increases, the normalized sum-rate difference 
becomes smaller and this hints that the approximation of the average sum-rate by the limiting sum-rate becomes 
more accurate. The difference is already about 1.3% and 0.5% for the MCP and CBf respectively for = 60, K = 
36. 

Figure 3 describes the applicability of the large system results into finite-size systems. We choose a reasonable 
system-size in practice, i.e., N = 10, A' = 6. Then, 250 channel realizations are generated. For each channel 
realization, with a fixed regularization parameter of the precoder, the optimal u, denoted by v*^, is computed. Then 
the resulting average sum-rate is compared to the average sum-rate that using u* from the large system analysis, 
i.e., (18) and (25), for different values of e. We can see that the normalized average sum-rate difference, i.e., 
E[|Rsum(^s)-ftnm(i^ )|] ^ j^^^ ^ pcak around 4% that can be considered as a reasonable value for the chosen 

system size. For MCP, it is less than 0.47%. To this end, our simulation results indicate that the large system 
results discussed earlier approximate the finite-system quite well. 

In the following, we present some numerical simulations that visualize the characteristics of the optimal u* 
(in the large system limit) and the corresponding limiting SINR for each cooperation scheme. We are primarily 
interested in their characteristics when the interfering channel gain e varies, as depicted in Figure 4. In general, we 
can see that for the same system parameters, the CBf scheme allocates more power to feed back the direct channel 
compared to the MCP. From Figure 4(a), we can see that for values of e ranging from up to a certain threshold 
(denoted by e^J = e^!^x ^^'^ for MCP and CBf respectively), the optimal is 1 : in other words, it is optimal in 
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this range for the BSs not to try to get information about the cross channels and to construct the precoder based 
on the direct channel information only. Effectively, the two schemes reduce to the SCP scheme when u* = 1: as 
a result, the same limiting SINR is achieved by both schemes. 

In Figure 4(b), we can observe a peculiar behavior of the limiting SINR of MCP, which we already highlighted 
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in the analysis of Section III-A. When y/l < ^r^, i.e. when v* = 1, the SINR is decreasing as e increases. After 
that the SINR is still decreasing until e reaches e"p and then increasing: this reflects the trade-off between and 
b^. Note that this initial decrease does not occur in the perfect CSI case where the SINR is strictly increasing in 
e for MCR Similar to the MCP case, we can see that the limiting SINR of CBf is decreasing in e when v* = \ 
(SCP). Moreover, it is still decreasing when both BSs perform CBf. 

IV. Quantized Feedback via Random Vector Quantization (RVQ) 

In this section, we will derive the approximations of the SINR for the MCP (2) and CBf (3) by analyzing them 
in the large system limit. We use these approximations to optimize the feedback bit allocation, and regularization 
parameter. This joint optimization problem can be split into two steps. First, we derive the optimal bit allocation, 
i.e., the optimal B^ = ^ and Be = jf. Plugging the optimal bit allocation back into the limiting SINR expression, 
we can then proceed to the second step where we obtain the optimal regularization parameter. At the end of the 
section, some comparisons of the limiting SINR and bit allocation values for the two schemes are illustrated. 



A. MCP 

-1 

l+9{/3,p) 

system limit, the SINR converges in probability to a deterministic quantity given by 



Theorem 3. Let pu,Q = (1 + e) ^a/N and g{(3,p) be the solution of g{P,p) = \p+ -, , . ) . In the large 



SINR^c,a = 7e9 /3, Pm.o) ^77— 7^5 ^> (26) 

7e + (1 +5(/3,p„,e))^ 

where 



7d(l+<:) 
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is defined as the effective SNR and 



d = — . (28) 

1 + e 



Proof: Refer to Appendix II-B. ■ 
Theorem 3 shows that the Umiting SINR is the same for all users in both cells. This is not surprising given 
the symmetry in their channel statistics and feedback mechanisms. Moreover, the only dependence of the limiting 
SINR on the bit allocation is via 7e, which itself is a function of d: d can be interpreted as a measure of the 
total quaUty of the channel estimates; In fact, given that and Be are constrained to sum up to Bt, d in (28) 
highlights a trade-off between increasing feedback bits for direct channel and cross channel. Comparing (13) and 

(26) , we can immediately recognize an identical structure between them. The effective SNR expressions (14) and 

(27) also share a similar construction, where (1 + e)d^ in (28) can be thought to be equivalent to uj^ + 

Now, we move tho the first step of the joint optimization i.e., determining the optimal bit allocation that 
maximizes (26). It is clear from (26) that Bd and Be contributes to the limiting SINR through d. It is easy to check 
that the Umiting SINR is an increasing and a convex function of d. Thus, maximizing SINR^pQ is equivalent to 
maximizing d, i.e. solving (cf. Eq. (28)). 



max yjl — Xd + e\ I -■ (29) 

Xa&[X,,l] V Xd 

where Xt = 2~^\ Bt = jf- and Xd = 2"^''. The solution of (29) is presented in the following theorem. 

Theorem 4. SINR^^^ is maximized by allocating Bd = — log2(a^^) ^^^^ to feed back the direct channel information, 
and Be = Bt — Bd to feed back the interfering channel information, where x*^ is the positive ( real) solution of the 
following quartic equation 

x\-Xtxl + {eXt)\xd-l) = Q. (30) 
Proof: The first derivative of the objective function over Xd is given by 



^5E[d] _ 1 / 1 I eXt 

dxd 2 I VI - Xd xl /l- ^ 



(1 + = ^ I + 1 (31) 



and lim.3;^^Xt = c«, lim2,.^_j.i = — oo. Moreover, the objective function is a concave function in Xd since 
,d'^E\d] ill, , ,/o 2 eXt leXtf XA^^^^i 



dx\ 2 1 2^ x\ 2 x\ V Xd, 

The stationary point, x^, is obtained by setting the derivative equal to and it is the non-negative solution of 

x\ - Xtx\^ (eXtf(xd-l) = 
Since the objective function is concave over Xd, then xjj gives the global optimum point. 
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Now, let us discuss how the optimal bit allocation vary with e. Since Xd = x*^ satisfies (30), then by taking the 
(implicit) derivative of (30) w.r.t. e, we have 

'4= .0, forX.<.^<l. 

This implies that as e increases, x'^ (B^) increases (decreases). This is consistent with the intuition that for higher 
e, more resources would be allocated to quantize the cross channel information. At one of the extremes, i.e., e = 0, 
x^ = Xt, or Bd = Bf. If e = 0, x^ = Xt, so that when there is no interference from the neighboring BS, all 
feedback bits are used to convey the direct channel states, as expected. At the other extreme, when e — )• oo, x^ — )• 1 
or Bd — )• 0. This can be shown by setting the derivative (31) equal to zero and we have 

1 _ Xt^/1 - Xd 
e " x^,,fl^' 

As e — > oo, the left hand side goes to zero and the stationarity is achieved by setting Xd = 1. 

It is also interesting to see how d, after optimal bit allocation, behaves as the cross channel gain varies. Let d* 
is d evaluated at Xd = x^. By taking we can show the following property. 

Proposition 1. For e < 1, d* is decreasing in e and increasing for e > 1. Consequently, d* is minimum at e = 1. 
As mentioned previously, x*, increases and consequently 1 — x*, decreases as e increases. On the other side. 



e^/l — Xt/x*^ is getting larger. So, from the calculation we can conclude that d* is mostly affected by ^1 — x^ 
for e < 1, while for the other values of e, the other term takes over. 

We now proceed to find the optimal Pmq that maximizes SINR^pg. The result is summarized below. 

Theorem 5. Let 7* be 7e evaluated at d = d*. The optimal Pm that maximizes SINR^„{d*) is 

Pl.a = ^- (32) 

le 

The corresponding limiting SINR is given by 

Proof: The equation (26) has the same structure as (13) and thus, (32) follows. ■ 
From Theorem 5, d* affects the regularization parameter and the limiting SINR via effective SNR 7*. The latter 

grows with d* (cf. (27)). Thus, /Jmq declines as the CSIT quaUty, d*, increases and this behavior is also observed 

for the cooperation schemes with the analog feedback. 

In Proposition 1 , we established how d* changes with e. We can show that 7* has a similar behavior but reaches 

its minimum at a different value of e due to the last term in the denominator in (27). For SINRj^'^ , it attains its 

minimum at e = eg, as described in the next proposition. 
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Proposition 2. Suppose that e = ej^ satisfies 



^ eXt [7d(l + e) + 1 + §] 

Then, SINRm'^q decreasing for e < and increasing for e > e". 

The characterization of SINR^'^q above reminds us a similar behavior of SINR^i^af after optimal power 
allocation. We can conclude that the limiting SINK of MCP under both feedback schemes has a common behavior 
as e varies. 

B. Coordinated Beamforming 

Theorem 6. Let p^g = a/N and Tq be the solution of the following cubic equation 

Pea + i+r^ + i+eFg 

Let (f)^ = I - 2'^", (j)c = I - 2'^", 5d = 2^^" and 5c = €2'^". In the large system limit, the SINE (3) for the 
quantized feedback via RVQ converges weakly to a deterministic quantity given by 

SW/?~ , = ^- ■ ^ (35) 



where 



(l+eFo)^ ^ (l+rQ)2 

Proof: See Appendix III-B ■ 
As in Theorem 3, Theorem 6 shows that that the limiting SINR is the same for all users. The quantization error 
variance of estimating the direct channel, Sd, affects both the signal strength (via (f)^) and the interference energy, in 
which it captures the effect of the intra-cell interference. 5c, on the other hand, only contributes to the interference 
term: It represents the quality of the cross channel and determines the strength of the inter-cell interference. Since 
Bt is fixed, increasing 3^, or equivalently reducing Be, will strengthen the desired signal and reduce the intra-cell 
interference: it does so, however, at the expense of strengthening the inter-cell interference. Thus, feedback bits' 
allocation is important in order to improve the performance of the system. 
To solve the joint optimization problem, it is useful to write (35) as follows 

SINR^jp = Gi- — , 

:^ + (1 - Xd){G2 - 1) + 6 (l - 1^ ) (G3 - 1) + 1 + e 



where Xd and Xt are defined as in the previous subsection and for brevity, we denote: Gi = —^'q (/^f^) ' ^2 = 
(1 + Fq)"^ and G3 = (1 + eP^)^^. The optimal bit allocation can be found by solving the following optimization 
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problem. 



max. SINR;?^f„. (36) 



The solution of (36) is summarized in the following theorem. 

Theorem 7. For a fixed Bt, the optimal bit allocation, in term of Xd = 2^^'', that maximizes SINR^^q is given by 



^ - 1-G3-Xt(2-G3) ~ 

€Xt{G3-l) + </e2Xt2(G3-l)2-eXt(^ + l+eG3)(G3-l) 

Xd = .^+1+60 — ~ ' otherwise. 

Id ^ 

Proof: Differentiating the objective function (36), we get 



(37) 



^""^ '4(i + (l-^^)(G2-l) + e(l-|f)(G3-l) + l + e 



2 



and the stationary is given by 



eXt(G3 - 1) + ^e^XfiGs - 1)^ - eXt{^^ + 1 + eG3)(G3 - 1) 



(38) 



Now let us consider the term Z = —x^{^ + 1 + eG^) + e(G3 — 1) (2X^2;^ — Xt) in the numerator. It can be verified 
that the sign of Z is the same as the sign of — dx™'^ ■ Thus, X^ = will be the unique positive solution of the 
quadratic equation Z = 0. 

It can be also checked that = — 2xrf(— + 1 + eG^) + eiG^ — l){2Xt) < and thus, Z is decreasing in Xd- 
Since at x^ = 1, Z < 0, we should never allocate x^ = 1. We will allocate x^ = if Z < at x^^ = Xt (this 
condition is satisfied whenever e < eth). ■ 

Unlike the MCP case where xj^ = Xt only when e = 0, in the CBf, it is optimal for a user to allocate all Bt to 
the direct channel when < e < eth. Note that xj = Xt does not imply that the cooperation breaks down or that 
both BSs perform single-cell processing. It is easy to check that eth increases when Bt or 7^ is decreased. This 
suggests that when the resource for the feedback bits is scarce or the received SNR is low then it is preferable 
for the user to allocate all the feedback bits to quantize the direct channel. So, in this situation, quantizing the 
cross channel does more harm to the performance the system. However, as e increases beyond eth, quantizing the 
cross channel will improve the SINR. We can show that x^, particularly Xd, is increasing in e. In doing that, we 
need to take the derivative of Xd over e. It is easy to show that Fq is decreasing in e. Then, it follows that G3 is 
decreasing in e. Using this fact, we can then show > 0. So, as in the case of MCP, this suggests that more 
resources are allocated to feedback the cross-channel when e increases. 

Once we have the optimal bit allocation, we can find the optimal pc,Q, as we did for the MCP. For that purpose. 
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we can rewrite (37) w.r.t pcg as follows 



Xt, PcQ > Ah 

Xd, otherwise. 



where for given Xt,e and 'yd, the threshold pth satisfies e = eth- So, we have SINR^fQ(Xd) for pcg < Pth and 
SINR^fQ(X() for other values of pco- 

Now, let us investigate the optimal pc.g when = Xd- By evaluating = 0, we can determine the 

stationary point, which is given by 



'- + Xd + e^ 
.Id Xd 

/5 



+ e(G2G'3 - GgG'^) 



X'a{G2 + eGsh, 



X> + (1 - Xd) ((1 - X,)G'2 + e{l-^) G's) 



where 7e 



7d 



l + e + e(G3-l) 1 



2Xt 



+ 4i I and G'n = #^ and 



dpc.Q ' 



We can show that the derivative is positive for pc.g G [OiPx^) and negative for pc.g G (Px^'^^)- Since 
SINR^fQ(Xd) is defined for pc,Q < pth, if P°Xa < Ah then SINR^fQ(Xrf) is increasing for pc,Q G [0,p°xj and 
decreasing for /Oc.q € [p^^^, Ah)- If > Ah then SINR^fQ(Xd) is increasing for />c.q S [0, Ah)- 



Then, we move to the case when = Xt. By setting 



dpc.Q 



0, the stationary point is then given by 



(G'2 + eG'3)iXt + l/7d + e) + e(l - Xt)iG2G'^ - G^G'^ 



(l-Xi)G'2 

We can also show that the derivative is positive for pc.g G [OiPxJ ^'^'^ negative for pc.Q € (p^^rjiOo). Since 
SINR^fQ(X() is defined for pc.g > Ah> if P°x^ > Pth then SINR^fQ(Xt) is increasing for pc.g G [Pth,PxJ ^'^'^ 
decreasing for pc.g G [Px^jOo)- If P°x^ < Pth then SINR^fQ(Xt) is decreasing for pc,g G [pth,oo). 

In what follows, by knowing the stationary point in both regions of p, we will investigate how to obtain the optimal 
Pc.g, denoted by p^o' for Pc,g £ [0,oo). By inspecting 0SINR^fQ(Xd)/9pc.g and 5SINR^fQ(X()/5pc,g we can see 
that that SINR^fQ(x^) is continuously differentiable for the region, pcg G [0,pth) and pcQ G [pth,oo), respectively. 
To show SINR^jq(x^) is continuously differentiable for pcg G [0, 00) we need to establish 9SINR^fQ(x^)/(9pc,g 
to be continuous at pc,g = pth, or equivalently 



^.^ aSINR-o(^d) 



Pc.g-s>P,h 



dp, 



aSINR-,(Xt) 
lim - 



aSINRS,Q(Xt 



'eg 



Pcq^pI 



dp 



eg 



dp 



c.g 



(39) 



Pc,g=P(h 



When Peg — )• P(jj, X^ — )• Xt and therefore the denominator of dS^NR^jQ{Xd) / dpc.Q and 9SINR^,Q(Xf)/9pc.Q are 
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equal. Let M{f) denote the numerator of /. As Xd — )• Xt, we have 

lim N (aSINRS,Q(Xd)/a/9c,Q) = [/3(G'2 + €G'^){lhd + 1 + e - (1 - X*)) + /3e(l - Xt){G2G'^ - G^G'^) 

PC,Q->P,h 

-Ah(l - Xt)G'2] Tq{1 - Xt) - lim X'^{pG2 + (ieG^ + phe 



( aSINR- 



\ dp^ 



C.Q 



lim x:i(/3G2 + /3eG3 + p)7e 



PC.Q=Pth / PC.Q->P,h 



where = dX^/dpcq- We should note that lim ^ - = — AeGg j_ , ^ ' 7^ 0. This shows that is not 

rC,Q r th H~ J- I ^ 

continuously differentiable over pc.Q- Now let us look at 



lim 7, = J- + l + e + e(G3-l)(-l + ^ 

Pc.o-->P,h 7d V 



1 

Xt 



- + 1 ) Xt + t{2Xt - 1) - eG3(Xt - 1) 
Id 







since as pco —5" Pth, from the (equivalent) condition e = eth, the term in the bracket becomes 0. This concludes 
(39) and therefore SINR^fQ(j;^) is continuously differentiable for pcQ G [0, 00). 

By using the property above and the facts that the SINR^fQ(Xd) and SINR^fQ(Xt) are quasi-concave (unimodal), 
we can determine the optimal and jointly as described in Algorithm 1. We can verify the steps 6-13 in 



Algorithm 1 Calculate p* and x 



9: 
10: 
11: 
12: 
13: 
14: 



Compute pth 
if Pth < then 

4 = Xt. 

* o 
PCQ = Px, 

else 

Compute p°-^^ 
if P°x^ > Pth then 
^*d = Xt 
Pc.Q = Px, 

else 

^d = Xd 

Pc.Q = Ph 

end if 
end if 



the algorithm by using the following arguments: If p°^^ > pth> then the derivate of SINR^fQ(Xt) is positive at 
Pc.Q = Ah because SINR^fQ(Xt) is quasi-concave. Since the SINR^jg(x^) is continuously differentiable, then 
the derivative of SINR^fQ(Xrf) is also positive when pcQ — )• pth- Since SINR^fQ(Xrf) is also quasi-concave, 
consequently SINR^fQ(Xd) is increasing for pcQ G [0, Pth)- This implies that p*^ = p\^. Similar types of arguments 
can be also used to verify that if p\^ < pth then = P°Xi- 
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5 10 15 20 25 30 35 40 45 50 55 60 



N 

Fig. 5. The total sum-rate difference for different system dimensions with /3 = 0.6, e = 0.5 7d = 10 dB and Bt = 4. 
C. Numerical Results 

The first two figures in this section are obtained by using a similar procedure to that followed in the analog 
feedback case. Figure 5 shows how well the limiting sum-rate (equivalently the limiting SINR) approximates 
the finite-size system sum-rate. The optimal regularization parameter and bit allocation are applied in computing 
the limiting and average sum-rates. As N grows, the normalized sum-rate difference become smaller. For = 
60, K = 36, it is arleady about 3.1% and 1.6% for MCP and CBf, respectively. Figure 6 shows the average 
sum-rate difference, with a fixed regularization parameter, between the system that uses B^^^ and Bd* to feed 
back the direct channel states. B^ pg denotes the optimal bit allocation of the finite-size system. For each channel 
realization, it is obtained by a grid search. With = 10, = 6, the maximum normalized average sum-rate 
difference reaches 0.22% for MCP. It is about four-times bigger for CBf, which is approximately 0.86%. Thus, 
from those simulations, similar to the analog feedback case, the conclusions we can reach for the limiting regime 
are actually useful for the finite system case. 

In the following, we present numerical simulations that show the behavior of the limiting SINR and optimal bit 
allocation for MCP and CBf as e varies. The optimal bit allocation is illustrated in Figure 7. As shown in Section 
IV, the optimal B^ for MCP is decreasing in e and B^ = Bt when e = 0. For CBf, B^ = Bt when e < 0.19, and 
after that decreases as e grows. Overall, for given e, B^ for CBf is larger than for MCP, implying the quality of 
the direct channel information is more important for CBf. 

In Figure 7, the optimal values for the regularization parameter and bit allocation are used. From that figure, it 
is obvious that SINR^fQ decreases as e increases. In the case of MCP, as predicted by the analysis, the limiting 



25 




Fig. 6. The (normalized) average sum-rate difference of the finite-size system by using the -B^j ^s ™d with N — 10, = 0.6, 7d — 10 
dB and = 4. 



SINR is decreasing until ej^^yg f« 0.72 and is increasing after that point. By comparing the limiting SINR for 
both cooperation schemes, it is also interesting to see that for some values of e, i.e., in the interval when CBf has 
B* = 0, the CBf slightly outperforms MCP. We should note that within the current scheme, when B* = 0, CBf 
and MCP are not the same as single-cell processing (SCP): under RVQ, there is still a quantization vector in the 
codebook that is used to represent the cross channel (although it is uncorrelated with the actual channel vector 
being quantized). 

Motivated by the above facts, we investigate whether SCP provides some advantages over MCP and CBf for 
some (low) values of e. In SCP, we use Bj^ jj = Bt bits {\/k,j) to quantize the direct channel. The cross channels 
in the precoder are represented by vectors with zero entries . By following the steps in deriving Theorem 3 and 6, 
we can show that the limiting SINR is given by 

l + f(l + <7(/3,Ps))2 



SINR^p„ = 7e5(/3,Ps 



7e + (l+5(/3,Ps))2 



where ps = N and 7e = 2-Bt'4-e+-L ■ follows that the optimal ps maximizing SINR^p^ is p* = — and the 
corresponding the limiting SINR is SINRs^^ = QiP^Pt)- 

From Figure 7, it is obvious that the SCP outperforms MCP and CBf for some values of e, that is, e < 0.13. 
Surprisingly, the CBf is still beaten by SCP until e 0.82. This means that the SCP still gives advantages over 
the CBf even in a quite strong interference regime with this level of feedback. 
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Fig. 7. (a) Optimal bit allocation vs e (b) Limiting SINR vs. e. Parameters: 7^=10 dB, Bt = 4. 



V. Analog vs. digital feedback 

In this section we will compare the performance of the analog and quantized feedback for each cooperation 
scheme. For the quantized feedback, we follow the approach in [14], [21], [23], [38] that translates feedback bits 
to symbols for a fair comparison with the analog feedback. In this regard, there are two approaches [14]: 

1) By assuming that the feedback channel is error free and transmitted at the uplink rate (even though this 
assumption could be unrealistic in practice), we can write 

Bt = ^ = 2Klog^{l + {l + ehu). (40) 

This approach is introduced in [21], [23]. (40) is obtained by assuming that each feedback bit is received 
by both base stations in different cells where the path-gains from a user to its own BS and other BS are 
different i.e, 1 and e respectively. We can think the feedback transmission from a user to both BSs as a 
Single-Input Multi-Output (SIMO) system. The BSs linearly combine the feedback signal from the user and 
the corresponding maximum SNR is (l + e)7u (see [39]). The pre-log factor 2kN for Bt in (40) presents the 
channel uses (symbols) for transmitting the feedback bits which are the same as those for the analog feedback. 
K follows the discussion in Section II-C. Our approach is different from the approach in [14] in which the 
user k in cell j sends the feedbacks only to its own BS j. In that case, (40) becomes Bt = 2Klog2 (1 + 7m)- 

2) Following [38], the second approach translates the feedback bits to symbols based on the modulation scheme 
used in the feedback transmission. In the analog feedback, the feedback takes 2kN channel uses per user. 
Let 77 be a conversion factor that links the bits and symbols and it depends on the modulation scheme. As 
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an example, for the binary phase shift keying (BPSK), rj = 1. Thus, we can write (see also [14]) 

r]Bt = 2kN. (41) 

We should note that using this approach, for a fixed k there is no link between Bt and 7^ as we can see in 
(40). 

Let us assume that k = 1. Thus, with the first approach, we have Xt = 2~^* = yj . The comparison 

of the limiting SINR based on the analog and quantized feedback for MCP and CBf can be seen in Figure 8(a). 
It shows that the quantized feedback beats the analog feedback in both MCP and CBf for e less than about 1. 
A similar situation still occurs for CBf even for e G [0,2]. The opposite happens for MCP when e is above 1.5. 
The comparison of the analog and quantized feedback with the second approach, also with k = 1, is illustrated in 
Figure 8(b). Similar to the previous, one can see that the quantized feedback outperforms the analog feedback if 
e is below a certain threshold. Otherwise, the analog feedback gives better performance. Those observations can 
be explained by verifying whether the feedback scheme that provides better CSIT will give a better performance. 
This is easier to check by looking at the MCP scheme because from our discussions in Section III and IV, its 
performance can be measured by the total CSIT quality, i.e., lOc + cod in the analog feedback and (1 + e)(P in the 
digital feedback. Plotting those over e, not shown here, will give the same behaviors for the MCP as we observed 
in Figure 8. Thus, from our simulations above, the CSIT quality of the quantized feedback is better than that of 
analog feedback when the cross channel gain is below a certain threshold. The above plots also confirm that more 
feedback resources will increase the system performance: for a fixed 7^ = dB, Bt in the left plot is larger than 
that in the right plot and hence gives a higher (limiting) SINR for the quantized feedback scheme. 
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Figure 9 depicts the limiting SINR of the analog and quantized feedback for different values of feedback 
rate. For the analog feedback, the values of feedback rate/bit is converted by using the previous approaches: 
K = 2 log ) ^i^d ^ = Bt/2 respectively. For MCP, we can see that initially the analog feedback scheme 

outperforms the quantized feedback in both plots. However, after a certain value (threshold) of Bt, the opposite 
happens. A similar observation also holds for the CBf scheme. The explanations for those phenomena follow the 
discussions for Figure 8. We should note that in generating the figures, the values for Bt are already determined. 
So, the limiting SINRs for the digital feedback are the same in both sub-figures. For the analog feedback, since 
K with the approach (41) is larger (with 7^ = dB) than that with the approach (40), then the training period in 
the former is longer and will result in a better CSIT. Thus, the limiting SINRs for the analog feedback in Figure 
9(b) are larger compared to those in 9(a). 

VI. Conclusion 

In this paper, we perform feedback optimization for the analog and quantized feedback schemes in a symmetric 
two-cell network with different levels of cooperation between base stations. In both cooperation schemes, it is 
shown that more resources, uplink transmit power in the case of analog feedback or feedback bits in the case of 
quantized feedback, are allocated to feeding back the interfering channel information as the interfering channel 
gain increases. Moreover, if the interfering channel gain is below a certain threshold, the conventional network 
with no cooperation between base stations is preferable. Our analysis also shows that the limiting SINR for MCP, 
in both analog and quantized feedback, improves in e if e is above certain threshold. This also implies that above 
that threshold the (total) quality of the channel at the base stations is also getting better Although our analysis 
is performed in the asymptotic regime, our numerical results hint to their validity in the finite-size system cases. 
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Future works could consider a more general channel model such as analog feedback through MAC channels. 
Furthermore, feedback reduction problems in which the users or groups of users have different path-loss gains can 
be interesting to explore. 



Appendix I 
Some Results in Random Matrix theory 

For the clarity in presentation, in this section we will list some results in random matrix theory that will be used 
to derive the large system results in this work. 

Lemma 1 ( [40, Lemma 1]). Let A be a deterministic N x N complex matrix with uniformly bounded spectral 
radius for all N. Let q = -^[qi,q2, - ■ ■ ,qN\^ where the Qi's are i.i.d with zero mean, unit variance and finite 
eight moment. Let r be a similar vector independent of q. Then, we have 

qAq^ - ^7>(A) ^ 0, and qAr^ ^ 0. 

Theorem 8 ( [41]). Let H be a \_cN\ x \_dN\ random matrix with independent entries [H]jj which are zero mean 
and variance E [|[H]jjp] = N^^Pij, such that Pjj are uniformly bounded from above. For each N, let 

VN{x,y) : [0,c] x [0, d] ^ M 

be the variance profile function given by 

i i+1 j J+1 

VN{x,y) = P,j, -<x<^, _<y<__. 

Suppose that vn{x, y) converges uniformly to a limiting bounded function v{x, y). Then, for each a, 6 G [0, c], a < b, 
and z G 

1 Y: [{UU^-ziy']^^^ u{x,z)dz 



i=\aN\ 

where u{x, z) satisfies 



u{x, z) = 



_^ _|_ v{x,y)dy 

Jo u{w,z)v{'w,y)d'w 

for every x £ [0, c]. The solution always exists and is unique in the class of functions u{x, z) > 0, analytic on 
z G and continuous on x £ [0, c]. Moreover, almost surely, the empirical eigenvalue distribution of HH^ 
converges weakly to a limiting distribution whose Stieltjes transform is given by u{x, z) dx 

In the theorem above x-axis and y-axis refer to the rows and columns of the matrix H, respectively. 
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Appendix II 

Large system results for the Network MIMO 



First, we will expand the SINR expression in (2). Let j = diag{(/)fc j i, i;^)^ j 2}- Based on (1) we can write 
hfc j = hf^j^^ J + hk j. Consequently, the SINR can be expressed as 



SINR 



k,j 



(42) 



-1 



It holds that ^H^H + alj = {Hj^jHkj + hj^jh^j + alj . By applying the matrix inversion lemma (MIL) 
to the RHS of the (previous) equation, we obtain 



H^H + al 



H^,H;.., + gi) h^,fi,, (nj^jU,, + gl) 



Let p = jj and Ofcj = ^ -^H^^^H^ ^ + pi j . Then, we can write (43) as -^Z^j, where 



Thus, (42) can be expressed as 



^k,j + Fk,j 



SINRfc 



1 + A 



c2 (S,,,+25ft[Dfc,,]+£;fc,,) + a2 



,2 ' 



where 



(43) 



(44) 



Bh 



Ek,j 



H 



^k,j = —hk,j^ljOk,jK-j 



hkjOkjhj^j 



N 
1 



1 



1 



Dk,j - —hkj^^jZkj l—H^jUkjj Zkjh^j 



N 



:^k,jZkj 



^k,j^k,j ) Zfc,jhf 



N 



H 



Note that for the analog feedback, ^k,j = I;VA;,j. In the following subsections, the large system limit of each 
term in the SINR (44), for the analog and quantized feedback cases, will be derived. For brevity in the following 
presentation, we also denote Qkj = Okj (jf^k^kj) Okj- 
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A. Proof of Theorem 1: Analog Feedback case 

In the analog feedback case, as mentioned in Section III, h^ j ~ CAA(0, WjjlAr) and hfc^jj ~ CAA(0, ^jjlAr) are 
independent. We can rewrite those vectors respectively as 

hfej = %kM,j and hk,j = dfcjD|_^., 

where j ~ CAA(0, l2Ar) and j ~ CAA(0, l2Ar) are independent. The diagonal matrices ,,■ and j are given 
by Gfcj = diag{a;jilAr,a;j2lAr} and D^j = diag{(5jilAr, (5j2lAr}, respectively. 

Since in the analysis below, we heavily use Theorem 8, let us define the asymptotic variance profile for the 
matrix ;^H^ which is a 2N x 2(5N complex random matrix. Following Theorem 8, in our case, we have x G [0, 2], 
y G [0, 2/3] and the asymptotic variance profile is given by 

ujd 0<x<l,0<2/</3 
ujc l<x<2,0<y</3 
ujc < X < l,/3 < y < 2/3 
ujd 1 < X < 2,/3 < y < 2/3. 
In what follows, we will derive the large system limit for the terms in (44). 

i i i 

1) Akj: It can rewritten as j^gkjG^j^'^ jOkjG'^ jgj^j. Applying Lemma 1 yields 



v{x,y) 



Akj - ^Tr ( Ck.,^l „.0 



0. 



The second term in the LHS can be written as 
applying Theorem 8, it converges in probability to 



N n \r\ ^ . v^2Af 



Ei=i ^ji\/4>k,j,i[Ok,j]ii + Ei=Af+i ^j2\/(t>k,j,2[Ok,j]i 



By 



u;ji\/(pk,j,i / u{x,-p) dx + ujj2V(pk,j,2 / uix,-p) dx, 



where for < x < 1, 



and for 1 < X < 2, 



u{x, —p) = Ui 



U{x, -p) = U2 



P + 



+ 



l+UlljJd + U2U)a l+U-i_UJa+U2U!d 



1 



„ I M"^c I PijJd 

" l + UiL0d+U2Ula l + UiU}a + U2iOd 

The solution of the equations above is ui = U2 = u where u is the positive solution of 

1 



P{uJd+OJa) ' 
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Let g{f3,p), be the solution of g{(3,p) = {p + 



1 



u 



/3 

l+9(/3,p) 



Then, we can express u in term of 5(/3, p) as 



We should note that Uji + Wj2 = + cjo Vj. Thus, 



- i^, ^ 0, with AT. - ^^^^-1 + 



Since in the current feedback scheme j = l2Ar then A^- = g{f3,p). Moreover, A'^j = A°° is the same for all 
users in both cells. 

2) Akj: This term is Akj with = l2Af- Thus, it follows that A^j — g{(3,p) — ^ 0. 

i i ^ 

3j Ffc j.- It can be rewritten as j^dk.j^l jOkj^l j^k j- Conditioning on H^j, it is obvious that j, gfc j and 
D^^Ofc jG^^ are independent of each other. By Lemma 1, it follows that j — )• 0. 
4) Dkj: Expanding D^j, we have 



1 + A 



1 + ^ 



I)(1)+Z)(2)_Z)(3)_I)(4). 

fe,J ' fe,J KJ fej 



(45) 



Following the arguments in 4), it can be checked that D^^j 0. Similarly, D^^j 0. Moreover, since F^j 
then it follows that d\^] and D^^] 0. Combining the results yields D^j 0. 



5j It can be rewritten as 

1 



= <]+ 



1 + A 



A* R*-^-* 



(46) 



(1 + Afcj)2 1 + Ak,j 

From Lemma 1, we can show B^^j — -^Tr (Gkj^kjQk.j) 0. It can be shown that the following equality holds 

(47) 



d 



(xx^ + ci)"^xx^(xx^ + ci)"^ = (xx^ + ci)"^ + c^(xx^ + ci)"^ 

Thus, we have Qkj = Ok,j + P^^k,j- By applying Theorem 8, we obtain 



[hN\ 



77 XI 



t.p. 



i=[aN\ 



^ d 

u{x, —p) dx + p— I u{x, —p) dx. 



Consequently, we can show B^p. - ^ 0, where 

^(l),oo _ l^jr'l^fcj,! + ^j24>k,j,2 



dp 



d 

g{/3,p) + p-^g{p,p) 



(48) 
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(2) „(2),oo 



Similarly, we can also show that ■ — Bf^ ■ 



(3) r(3),oo 



and B):' - B):'. 



0, where 



<]'°°=5(/3,p) + m(/3,/^) 



dp" 



d 



For the analog feedback case, i.e., j = '^2N, it follows that B 



B 



{2),oo 



j(3),oo 



. Therefore, S^^j 



converges in probability to 



g{P,p) + p-QzgW,p) 



6) Ekj: Expanding this term gives 



Ek,j = —hkjQkjh^j - 23f? 



^k,j ^k,j ^ ^k,j ■ 



1 + ^fcj 



+ 



Fk,j\ -w^k,jClk,j^k,j 
(1 + 



(49) 



By using the previous results, we can show that E^^'^- and E-^ 



,(3) if- 



0. From Lemma 1, it follows that 



E^^j — ;^Tr (Dyt jQfc j) 0. Following the steps in obtaining B)^'^'^ , it is easy to show that E)^'- —Ej^ ■ 
where 



,{l),oo 



(1) jj,(l),oo i.p. 



0, 



E: 



k,j 



Thus, we have Ekj - e'^]'"" ^ 0. 



(l),oo _ ^k,j,l + SkJ,2 



_d 

9W,p) +p-^g{l3,p) 



7) c^: The denominator of can be written as follows 



i^Trf fi-H^H + pI^ '^H^H 



N' 



N 



N 



where i^Qffg is the empirical eigenvalue distribution of H^H. From Theorem 8, ^gng converges almost surely 
to a limiting distribution G* whose Stieltjes transform m{z) = J^^ dG*{X) = u{x,z) dx. Therefore, 

A 



1 9 1 

X + p dpX + p ^ 

H> m{-p) + p-^m{-p) 
d 

/ u{x, -p) + p—u{x, -p) dx. 
Jo dp 



(50) 



Previously, we have shown that u{x, —p) = 2{ujd + ^c) ^g{f3, p) with p = 2p{ijjd + (^c) ^- Hence, (50) is equal 



to 2{ujd + ojc) ( (7(/3, p) + p-^g{(3, p) ) and the large system result for c is given by 



\{u:d + uJc)Pt 

g{P.p) + Pmg{l^^p) 



0. 
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The large system results in 1) - 3) and 7) show that the signal strength, i.e, the numerator of (44), converges to 

iud + u}c)Pg\/3,p) ^^^^ 



where we already substitute P = \Pt. Similarly, it follows that the interference (energy) converges to 

Ud + uc- {6d + 6c) (1 + P)f) ■ 



(52) 



{l + g{P,p)f 

To simplify the expression for the limiting SINR^ j, we need the following result (see [37]) 

9{P, P) + PTp9{P, P) = n , , to 

dp /3 + p{l + g{f3,p)y 

By combining the large system results and also denoting /j^ af = P> we can express the limiting SINR as in (13). 
This completes the proof. 



B. Proof of Theorem 3: Quantized feedback (via RVQ) case 

Since the 4'k,j,i is a random variable, while it is not in the analog feedback, the derivation for the limiting SINR 
is quite different from the previous subsection. For a given 2N x 2N matrix X, we can partition it as follows 



^k,j ^k,j 
-v-21 v22 



21 



Where X^i. = [X^,,-],^, Z = 1, • • • , iV, m = 1, • • • , iV, X^^. = [Xfc,,]^™, 1 = 1,... ,N,m = N+1, • • • , 2Af , X 
[Xfcj],„,Z = A^ + !,••• ,2iV,m = !,••• ,A^, and X22. = [Xfc,,-],„, / = + 1, • • • , 2iV, m = iV + 1, • • • ,2N. 

In the following, we derive the large system limit for each term in the SINR. 

1) Akj: We can write A^j as 

Ak,j = (^(f>^ j^^h.k,j,iOl]jh^j i + 4>^ j^^hkj^iOl^jh^j^2 + ^k,j,2^k,j,20l]jhj^j^i + (f>^j2^k,j,20l^jhj^j^2) ■ (53) 

Since hkj^i and h.kj,2 are independent then the second and third terms converge almost surely to 0. For the 
first term 



1 



^K.iO^!,h-,-^Tr(0-.) 







or equivalently. 



1 C rjii v^^ 



N 



N 



i=l 



By using the Girko's result (Theorem 8), we have 



N 



N 



k,j\ii 



l.p. 



i=l 



1 + e 



giP,p), 
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where p = By using the same techniques as in [16], we can show that 



1 - t2 



\J\ — 2-^'' j = i 
\/l — 2~^= otherwise. 



We should also note that the convergence in mean square sense implies the convergence in probability. By doing 
the same steps as for the last term of (53), we have 



A 



k,j 



1 + e 



2) Akj: This term is A^j with j = l2Ar. Hence, — g{P,p) 0. 

3) Fkj: We can expand the term as follows 



Fk,j - (hk,j,iOljhkj^i + hkj,iOl^jh^j 2 + ^k,j,20ljh^j^i + h^j- 20|jh^j- 2-) (54) 

By using the same steps as in deriving (65), it follows that each term on the RHS of the equation above converges 
in probability to 0. Hence, Fkj 0. 

4) Bk j.- By following the representation (46) for Bf^j , we can express the first term on the RHS as 



^kj = ]Y (^0fcj,lhfcJ,lQfejhf j_i + (j)k,j,lhk,j,lQijhj!j^2 + 'Pk,j,2^k,j,2QkJ^k,j,l + 4>k,j,2iik,j,2QkJ^k,j,2 

The second term and the third term of the equation above converge (almost surely) to 0. For the first term, 

lh.,,Q-h-,-^Tr(Q-0^ 
From (47), we have Q^ ., = O^j + p-^Okj- Hence, we can show 



1 



a2 uH 



»21 tH 



\22 tH 



N 



I. p. 



N 



i=l 



d 

u{x, —p) dx + p— I u{x, —p) dx 



dp 



1 + e 



g{f3,p)+p-g^g{f3,p) 



(1) 



By doing the same steps for the last term of Bj^ ■ , it follows that 



(1) 1-2-^^ + 6(1-2- 



B 



Similarly, we can also show that 



Bf^^-g{P,p)+p^_g{P,p)^Q, 



1 + e 



9W,p)+p-Qzg{/3,p) 



i.p. 



a/i - 2--Sd + eVl - 2-^-= 



d 



I. p. 



0. 
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Combining the results together, we obtain 

/ 1 _ 2-^'' + e(l - 2-^" 



k,j 



1 + e 



d?{2 + g{l3,p))g{l3,p) 
(l+5(/3,p))2 



_d 



l.p. 



where 



1 + e 



5) Dkj: We can expand D^ j as expressed in (45). Using the previous results, it follows immediately that D 

(3) 

and I?^ J converge to 0. We can expand the first term as follows 



(2) 



^kj - 1^ [^k,j,l^k,j,lQkj'^k,j,l + '^fcj- lhfcJ,lQfcjhfj- 2 + 4>lj^2^k,j,2QkJ^k,j,l + <^fcj,2hfc,j,2Qfcjhf J 2 



N 



»11 ZH 



,12 ZH 



21 ZH 



\22 ZH 



(1) 7^ 
k,j 



Again, by following the same steps as in deriving (65), each term converges in probability to and hence D 
0. Similarly, ^ 0. Therefore, D^j ^ 0. 

(2) (3) 

6) Ekj: The expansion of E/^j follows (49). By applying the previous results, £'^j and converge in 
probabiUty to 0. E^j can rewritten as 

-^fcj = ^ (^k,j,lQk]j^k,j,l + hA:,i,lQfejh^j,2 + ^k,j,2Q,l]j^k,j,l + ^k ,j ,20, l^j^^kj ,2^ ■ 

Since i and j 2 are independent then the second and third term converge to 0. The first term in the equation 
above can be written as 

H 



1 



^ V, nil r// _ '^fc,j,ill^fcJ.ill^ 

n-2 111, l|2 



(vfc,i,ih|^j.i)hfe 



(vfcj-ih|^.i)hfcj-i 



k,j, 



VfcJ,lQfcj-v|[j- + 



J_llli l|4 

JV2 ll"A;,j,l|| 



-23f? 



J_I|C 112 



Since Vk,j,i and 1 are independent, the second and third term in the bracket converge to 0. It can be shown 
that 



1 oil 1 

j;^^k,j,iOk,j^k,j,i - Y 



+ e 



l.p. 



0. 



The large system limit for the last term of Ek,j can be done in the same way. Thus, 



Ekj 



+ e2- 



1 + e 



_d 

9{l3,p) + p-gz9W,p) 



l.p. 



0. 



7) c^: We can show that 



0. 
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Since implies then also converges weakly to the same quantity as above. 
Combining the results, it follows that the signal strength and the interference converge to 



and 

respectively. Therefore, (26) follows immediately with Pm.q = P- 



Appendix III 

Large System Results for the Coordinated Beamforming 
For brevity in the proofs, we define the following (see also [17]) 

2 K 



\ m=l 1=1 J 

V («,»n)^(fej) / 



-1 



^kj,k'j',j 



{l,m)y^{k,j),{k',j') 

From the definitions above, we can write the numerator of the SINR^j (3) excluding c^, as follows 



\^k,j,j^kj\ 



\l (t>k,j, 

'^^k,j,j-t^kjilk,j,j 



+ 



1 

N' 



^k,j,j-^kj^k,j,j 



+ 2K 



1 

iV2 



(hfc, j ,i Afcj hj^jj ) (hfe J J Akj h^j^ 



5, 



(1) 

kj 



(55) 



(56) 



In denominator, let us consider the term |hfc j jvWfc/jvp which can be expanded as follows 

\K,j,j'Wk'j'f = j-^4>kjj'Kjj'A.k'j'hk\j'j''^k'j\j'^k'j''^k^^ + j-^'^kjj'-^k'j'iik'j'j'K'j'j'Ak'j'h^jj, (57) 



23? 



"fcjj'^fc'i'"fc'j',j'"fc'.i'j'^fc'i'"fe,jj' 



r(2) 



(3) 



(58) 
(59) 



jY kj,k'j' ^ kj,k'j' kj,k'j' jy 

Now, we are going to derive the large system limit for -^Tr (A^) since it will be used frequently in this section 
Let Hj = [hi,ij • • • hx.ij i^i,2,j ■ ■ ■ iiK,2,j]'^ and hk,ij ~ CM{0,u}ijl). Then, Aj = (p + -^HfHj) and 



N ^ J X + p 
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where i^pffp is the empirical eigenvalue distribution of H^^Hj. From Theorem 8, this distribution converges almost 
surely to a limiting distribution F whose Stieltjes transform 'mf:{z). It can be shown that 



1 as 

— Tr (Aj) mF(-/>) = J u{x, -p) dx. 



where 



uix,-p) = u{-p) 



for < X < 1. Let T = u{-p), then ^Tr (Aj) ^ T. 



A. Analog Feedback 

Based on the channel model (1), we have (l)k,j.j = (Afcj.j' = 1- The definitions for other terms such as uj, and 5, 
can be seen in Section II-C. Now, let us first derive the large system limit for the numerator of the SINR^j. We 
start with the term From Lemma 1 and by applying [42, Lemma 5.1], we can show that 



max 

j=l,2,k<K 



0. 



By applying rank-1 perturbation lemma (see e.g., [17, Lemma 3], [43, Lemma 14.3]), we have 



max 

j=l,2,k<K 



ail) _ ^^XrfA- 







where ^Tr(Aj) ^ T. 

Since h.kjj, Akj and hkjj are independent then it follows that 



max 

j=l,2,k<K 



Consequently, 



In summary. 



max 
j=l,2,k<K 



o{2) 



^k,j,j-^kji^kJJ 



and max 

j=l,2,k<K 



S', 



(3) 
kj 



max llhfcjjWfcj l - w^F 



2t2\ a.. 



0. 



j=l,2,k<K 

Now, let us move in analyzing the interference term. By using the matrix inversion lemma, we can rewrite 



hj,k'r as 



j(i) 

kj,k'j' 



N^k,j,j'^k,j,j'-^k'j',kj,j'iik' ,j' ,j'^k',j',j'-^k'j',kj,j'i^^jj' 



1 + ^hkjj'Ak'j'^kjj'^k,j,j'^ 
By applying Lemma 1, [42, Lemma 5.1] and rank-1 perturbation lemma (RIPL) twice, we can show 



max 

jJ'=l,2,k,k'<K,{kJ)^{k',j') 



1 

N' 



"^k,j,j'-^k'j',kj,j'ii-k,j,j' ^ Tr(Aj') 



LOii 



jj' ' 
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Similarly, 



max 

j,j'=l,2,k,k'<K,{k,j)^(k',j') 



jy'^'-k,jj'-^k'j\kjj'llk\j\j''^'-k\j'j'-^k'j\kjj'll-k,j,j' jy 



Since ^Tr (A^v) r and ^Tr (A 



1^, we have 



max 

j,j'=l,2,k,k'<K,(k,j)j^{k',j') 



kj,k'j' 



ojjj' dT 



0. 



By following the same steps, we obtain 



max 

j,j'=l,2,k,k'<K,{k,j)f^{k',j') 



(2) 

kj,k'j' 



d~p) 



0. 



and 



max 

j,j'=l,2,k,k'<K,{k,j)ft{k',j') 



j(3) 
kj,k'j' 



0. 



Combining the results, we have the large system limit for X in (59) 



max 

j,j'=l,2,k,k'<K,(k,j)jtik',j') 



I-UJd 



:i + ^,yr)2 



A 1^ 

' dp 



Using (60), the large system result for the interference term can be written as follows 

K K 



1=1 



UJr 



+ Sd + S, 



dp ' 



Now, we just need to derive the large system limit for cj = P ^X^a^i ll'^fcjlP^ > where we can express ||wfcj 



^hkjjAljh^jj. We can show that 



max 

j=l,2,k<K 



■^"^kjjAljh^jj - ^Tr (A^) 



Thus, 



where we can show that 



2 a..s 



P 



-PoJd'^ 



dp 



dp ^ 



r 



To sum up, from the analyses above, we can express the limiting signal energy as 



^Pu:dT[p + 



l3uJd 



(l+a;er)2 + (l+a;,r) 
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and the limiting interference energy as 

Finally, the limiting SINR can be expressed as (21), with = T and Pc af = P- 



B. Proof of Theorem 6: Quantized Feedback via RVQ 

In the derivation of the large system limit SINR in this section, we use some of the results presented in the 
previous section. Here, we have LOjj = lo^ = I and lo^j = lOc = e- From (1), we have 4>k,j,i = ^ ~ '^k j i- 

First, let us consider the numerator of the SINR. By using the result from previous section, we have 



max 

j=l,2,k<K 

where ■^Tr(Aj) T and T is the solution of 

r = 

As stated in [16], we have. 



^kj ^ l^J 



P+TTT + 



i+er 



h,, A 1 - 2-^^ 



(64) 



Since almost sure convergence and convergence in mean square imply the convergence in probability then 



^k,j,j\Sp,j 



(1)|2 



(1-2- 



0. 



By using (8), the term hkjjAkjhj^jj in 5*^^.^ can be rewritten as 



N 



H i7f^k,j,j^k,j,j)^k,j,j-Akj ^kj,j 



Since "Vk.jj and i^kjj independent, then 



max 

j=l,2,k<K 

It can also be shown that 



N 



0, and max 

j=l,2,k<K 



— VkjjAkjh^jj 



0. 



max 

j=l,2,k<K 



N' 



\h 



k,j,j I 



0, and max 

j=l,2,k<K 



N 



0. 



Hence, 



—hkjjAkjhkjj ^ 0, 



(65) 
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and thus, 



5g) ^ 0, and 4J ^ 0. 



Putting the results together, we have the following for the numerator 



Now, let us consider the interference terms. By using the same steps as in the previous section, we can show the 
foUowings 



max 

j,j'=l,2,k,k'<K,{k,j)jt{k',j') 



(1) _/ (i-2-^..'Hi'5r 

kj,k'j' \ 



max 

j,j'=l,2,k,k'<K,{k,j)ft{k',j') 



(l+^,yr)2 dp 

dp 



I. p. 



0, 



^kj,k'j' ^JJ''^ 



i.p. 



and 



max -'i./i-'o' — > 0, 

j,j'=l,2,k,k'<K,ik,j)^{k',j') "^''^ ^ 



where Bjji = when j = j' and otherwise Bjj' = Be. 
Combining the results, we have 



max 

j,j'=l,2,k,k'<K,{k,j)ft{k',j') 



(l+a;,yr)2 



dp 



i.p. 



Using (66), the large system result for the interference term can be written as follows 

K K 



(k'd'mkj) 



l=l,l^k 



1=1 



(i + r)2 + (i + er)2 i dp- 



dr 



By using the result from the previous results straightforwardly, we have ziwsl- Putting all the large 

system results for each term, we can show that the limiting signal strength is 



(66) 



(67) 



and the limiting interference energy becomes 

P 



+ 



+ Sd + S, 



(68) 



(i + r)2 (i + er)2 

Let pc Q = p and Fq = F. Then, we can obtain the limiting SINR given by (35) from (67) and (68) straightforwardly. 
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